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Abstract 

Diffuse ultrahigh energy 7-radiation can arise from a variety of astrophysical sources, including 
the interaction of 10 20 eV cosmic rays with the 3 K microwave background radiation or the collapse 
of topological defects created in the early Universe. We describe a sensitive search for diffuse 7- 
rays at ultrahigh energies using the CASA-MIA experiment. An isotropic flux of radiation is not 
detected, and we place stringent upper limits on the fraction of the 7-ray component relative to 
cosmic rays (< 10 -4 ) at energies from 5.7 x 10 14 eV to 5.5 x 10 16 eV. This result represents the 
first comprehensive constraint on the 7-ray flux at these energies. 
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We have known about the existence of ultrahigh energy (E> 10 14 eV) cosmic rays for more 
than 50 years 0], and yet their origin remains mysterious. Similarly, the possibility of a 7- 
ray component of the cosmic radiation remains an important open question in astrophysics 0|. 
Direct detection of ultrahigh energy photons by satellite and balloon-borne experiments is difficult 
because of the low fluxes. Fortunately, detection can be accomplished by large area (> 10 5 m 2 ) 
ground-based detectors using the air shower technique. This paper reports a sensitive search for a 
high energy diffuse 7-ray component of the cosmic rays made by the Chicago Air Shower Array- 
Michigan Muon Array (CASA-MIA) experiment. 

There are a variety of possible sources for a diffuse radiation at ultrahigh energies. Conven- 
tional astrophysical scenarios include: (a) unresolved point sources, such as active galactic nuclei 
(AGN) which may be the source of extremely high energy (E> 10 18 eV) cosmic rays ||, (b) decays 
of neutral pions produced by the collisions of cosmic rays with interstellar gas and dust ||], and 
(c) electromagnetic cascades resulting from the interactions of extremely high energy cosmic rays 
with the cosmic microwave background radiation (CMBR) ||. In scenario (b), the diffuse 7-ray 
flux is expected to be concentrated in the galactic plane. We have already reported a separate 
search for a diffuse galactic 7-ray component ||. 

In addition to these conventional sources, there might be new astrophysical or cosmological 
phenomena giving rise to a diffuse radiation at ultrahigh energies. For example, there is consider- 
able speculation that the highest energy (E> 10 20 eV) cosmic rays may result from the collapse or 
annihilation of topological defects formed at grand unified energy scales in the early Universe Q. 
With such scenarios, the production of cosmic rays is accompanied by electromagnetic cascades 
leading to copious 7-ray production. 

In any scenario in which 7-rays are produced at cosmological distances, we expect a substantial 
attenuation of the signal as a result of pair-production with intergalactic radiation fields. 7-rays 
near 10 15 eV should interact with the CMBR with an attenuation length ~ 20 kpc || . Nevertheless, 
cascading from higher energies in the presence of a weak intergalactic magnetic field can propagate 
7-rays over significant distances ||]. At the lower energies accessible by ground-based experiments 
(10 12 — 10 14 eV), absorption by infrared (IR) radiation is also expected, and may be significant, 
although the exact amount of absorption is still very uncertain ||. 

Ultrahigh energy photons and cosmic rays arrive at Earth, interact in the atmosphere, and cre- 
ate extensive air showers (EAS). EAS detectors sample the secondary shower particles (primarily 
electrons and muons). In showers initiated by cosmic-ray nuclei, muons are generated predom- 
inantly from the decays of charged mesons produced in hadronic interactions. 7-ray showers, 
however, are largely electromagnetic in nature since the cross section for meson photoproduc- 
tion is much smaller than that for pair production, by a factor of approximately 10~ 3 [10|. As 
a result, we expect a much smaller fraction of secondary hadrons in 7-ray showers relative to 



cosmic-ray showers, and correspondingly far fewer muons. Simulations based on the MOCCA [11 



and CORSIKA [12] programs indicate that 7-ray showers contain, on average, 3-4% of the num- 
ber of muons as showers initiated by protons of the same energy. These results, which agree with 
other calculations [13], are based on conventional electromagnetic and photoproduction couplings. 
The possibility of an anomalously high photoproduction cross section has been experimentally 
ruled out at HERA up equivalent laboratory energy of ~ 2 x 10 13 eV [ 10 1 . The simulation cross 
sections are in agreement with these findings, and are assumed to continue smoothly at higher 
energies. The technique for distinguishing 7-rays from cosmic rays by identifying muon-poor EAS 
was suggested more than thirty years ago [fbf| , and its main limitation results from downward 
fluctuations in the muon content of hadronic showers. In all experiments prior to CASA-MIA, the 
muon detectors have not been sufficient in size to provide a clear statistical separation between 
the 7-ray and cosmic-ray regimes. 
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Detections of ultrahigh energy 7-rays using the muon-poor air shower technique were reported 
in the 1960's by a Polish-French collaboration 15] and by the BASJE collaboration [16|. Later 
detections were reported by the Tien Shan [17] and Yakutsk [18] experiments. The first significant 
upper limit on the diffuse photon flux at ultrahigh energies was reported by the Utah-Michigan 
experiment using a portion of the eventual CASA-MIA detector [fl9|] . Upper limits have been 
reported by the HEGRA experiment, based on the study of the lateral distribution of Cherenkov 
light in EAS |20|] , and by the EAS-TOP detector based on muon-poor air showers |2l[| . The 
latter result used an analysis procedure in which the limit value was apparently optimized by an 
arbitrary cut on shower size. An upper limit based on the electromagnetic component of EAS 
recorded in emulsion chambers has been claimed p^ ], but this result is controversial p0| . 

The CASA-MIA experiment is located in Dugway, Utah, USA (40.2° N, 112.8° W) and consists 
of a surface array of 1089 scintillation detector stations enclosing an area of 2.3 x 10 5 m 2 , and a 
muon array of 1024 scintillation counters having an active of 2,500 m 2 |||. The muon array 
is more than ten times larger than any other EAS muon detector built for gamma-ray astronomy. 
For this analysis, we use data taken between January 1, 1992 and January 7, 1996, when the full 
experiment was operational. 

For each event recorded by CASA-MIA, we estimate the shower size from the number of 
surface detector stations struck, the shower core from the location of maximum particle density, 
and the shower direction from detector timing information |23| . The number of detected muons is 
determined from those muon counters hit within a narrow window (width ~ 100 nsec) around the 
expected arrival time. On average, we expect 0.65 muons per event from accidental muon hits. 

We have developed a set of data quality cuts, described elsewhere [24], to remove periods of 
time in which detector problems could potentially bias the reconstruction procedures. We require 
individual events to have valid information from the underground muon array, a shower core 
greater than 15 m from the physical edge of the surface array, and a reconstructed zenith angle 
less than 60°. We remove events in which a mismatch was likely to have occurred between the 
information derived from the surface array and that derived from the muon array. The fraction of 
these mismatched events is 0.034% of the total event sample, and is independent of shower size. 
After all cuts, the final data sample consists 1.370 x 10 9 events, which is larger by more than an 
order of magnitude than those used in other recent searches [19, [2(], 21 1. 

The median 7-ray energy for EAS that trigger CASA-MIA is ~ 10 14 eV. At this energy, the 
mean number of detected muons is ~ 7.5, and 4.1% of the events have zero detected muons. 
These zero muon events result from downward statistical fluctuations in the number of detected 
muons in ordinary cosmic-ray EAS. To greatly improve the background rejection, we consider data 
samples at higher energies in which the average number of detected muons is larger. We select 
five non-independent data samples by requiring the minimum number of stations to be greater 
than 50, 100, 250, 500, and 700. The total number of events, A^tot, and the average number of 
detected muons for each sample is shown in Table [l|. 

We use a simulation to reproduce the observed muon number distribution and to estimate the 
median cosmic ray energy of each data sample, as shown in Table |]. A library of artificial EAS are 
generated by the CORSIKA simulation code |l2|] using the known cosmic-ray energy spectrum (see 
discussion in |25|]) and a chemical composition that is independent of energy and consistent with 
spacecraft measurements [26|. The experimental detection efficiency is determined by observing 
the deviation of the shower size spectrum from a power law. Using the simulated muon size and 
assuming a Greisen lateral distribution function [27|, the number of detected muons is determined 
by a detector simulation which accounts for all known parameters of the muon detector (accidental 
muon hits, detection efficiency, etc.). We vary the median cosmic-ray energy in the simulation to 
obtain the best match between the simulation and data distributions of the number of detected 
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Table 1: Data samples used. The minimum number of stations required is listed in the first 
column. The total number of events, iVtotj the average number of muons, < Nn >, and the 
median comsic ray energy, E cr are given in the second, third, and fourth columns, respectively. 



Sample 


Ntot 


<N fl > 


E cr (TeV) 


> 50 stations 


6.9090 x 10 Y 


40.7 


575 


> 100 stations 


1.6042 x 10 7 


78.5 


1,350 


> 250 stations 


1.1863 x 10 6 


210.6 


5,000 


> 500 stations 


71,534 


470.1 


22,000 


> 700 stations 


11,572 


641.3 


55,000 



muons. The simulation reproduces the main features of the data for all samples. There is good 
agreement between the cosmic-ray energy values estimated from the simulation and those derived 
using the constant intensity method p5(| . 

We use the simulation to estimate the number of muons expected for showers initiated by 7-ray 
primaries. Fig. [l| shows the expected distributions in comparison with the data. For each sample, 
we define a signal region of the muon distribution in which the bulk of the 7-ray signal would be 
contained, and in which the known backgrounds are largely excluded. For example, in the > 250, 
> 500, and > 700 station samples, we exclude the bins at low muon content where we know there 
is potential background from mismatched events. No significant excesses consistent with a 7-ray 
signal are seen in any sample. Therefore, we assume that all events within the signal regions are 
background and set upper limits on the 7-ray flux and on the 7-ray fraction of the cosmic rays. 
We estimate, iVgo, the 90% C.L. upper limit on the number of detected events, using standard 
statistical methods [28], and use the simulation to evaluate the efficiency for 7-ray detection, e 7 , 



when the muon cut is applied. The upper limit on the fraction of the 7-ray integral flux relative 
to the cosmic-ray integral flux, J 7 /J cr , is given by: 

I cr e-yiVt t y Ey J 

where E cr is the median cosmic-ray energy, E~, is the median 7-ray energy, and a is the integral 
cosmic-ray spectral index (a = -1.7, E< 3 x 10 15 eV and a = -2.0, E> 3 x 10 15 eV). The factor 
involving energies in Eq. [l] accounts for the fact that 7-ray primaries produce larger showers than 
cosmic-ray primaries of the same energy, or conversely, that 7-ray showers are detected with the 
same efficiency as cosmic-ray showers of higher energy. Therefore, to report a flux ratio at a fixed 
cosmic-ray energy, we correct the 7-ray flux at its energy to the higher cosmic-ray energy under 
the assumption that both species have the same spectral index. Table ^ lists the values of iVgo, 
e 7 , and I^/I C r for each data sample. We use measurements of I cr to determine upper limits on 
the diffuse 7-ray flux at fixed 7-ray energies, using the values of iVgo, iVtot, and e 7 . The flux limits 
are given in Table [21 

Fig. |2| shows a compilation of measurements on the 7-ray fraction as a function of energy, 
including our work. The limits presented here represent the first stringent results spanning a 
wide range of energies, including the first results above 2 x 10 15 eV. The possibility of a diffuse 
7-ray flux at the 2 x 10 -4 to 4 x 10~ 3 level, as suggested by earlier results|l5], 16, 17, p], is ruled 



out. Our limits also constrain, to some extent, models @ which invoke topological defects as the 
cause of cosmic radiation at the highest energies. The limits disfavor models which have decaying 
topological defects at relatively nearby distances (i.e. < 100 Mpc), but not those which arrange 
them in cosmologically uniform distributions B. 
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Table 2: Results of the search for diffuse ultrahigh energy 7-rays. The median cosmic-ray energy, 
E cr , and the median gamma-ray energy, E 7 , are given in the first and fifth columns, respectively, 
in units of TeV. The quantities Ng and e 7 are defined in the text. I 1 jl cr is the 90% C.L. upper 
limit on the integral 7-ray fraction and 7 7 is the 90% C.L. upper limit on the integral 7-ray flux, 
in units of photons cm' 2 sec" 1 sr _1 . 



TP 


^90 




I 7 // cr (xl0- 5 ) 


-Ely 


/ 7 




575 


7273.0 


0.397 


< 10.0 


330 


< 1.0 x 10- 


-13 


1,350 


1718.0 


0.641 


< 6.5 


775 


< 2.6 x 10" 


-14 


5,000 


89.5 


0.663 


< 4.0 


2,450 


< 2.1 x 10" 


-15 


22,000 


2.3 


0.736 


< 1.5 


13,000 


< 5.4 x 10" 


-17 


55,000 


2.3 


0.891 


< 8.0 


33,000 


< 5.3 x 10" 


-17 



In summary, based on a very large sample of air shower events whose muon content has 
been measured by a large muon detector, we place stringent upper limits (< 10~ 4 ) on the 7- 
ray fraction of the cosmic rays at energies between 5.7 x 10 14 and 5.5 x 10 16 . These limits 
rule out the possibility of a significant component of the cosmic-ray flux whose interactions are 
predominantly electromagnetic in nature (e.g. 7-rays, electrons, positrons, etc.). In addition, by 
using measurements of the cosmic-ray flux, we place upper limits on the isotropic flux of 7-rays 
at these energies, which is of fundamental astrophysical interest. 
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Figure 1: Histograms of the number of detected muons for four data samples: (a) > 100 stations, 
(b) > 250 stations, (c) > 500 stations, and (d) > 700 stations. The data are represented by 
the points with error bars. The error bars account for the statistical uncertainty in the numbers 
of events and the systematic uncertainty in the subtraction of mismatched events. The 7-ray 
simulation results are represented by the solid lines, normalized to 10~ 3 of the cosmic-ray flux. 
The hatched regions correspond to the 7-ray signal regions as described in the text. 
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Figure 2: Measurements of the fraction of 7-rays relative to cosmic rays at ultrahigh energies. 
The hatched region indicates the range of 7-ray detections as reported in the literature prior to 
1985 [15,16]. The results of the the Tien-Shan [17] and Yakutsk [18] experiments are shown by 
the point with error bars and the asterix, respectively. The points with arrows represent upper 
limits from the Utah-Michigan [19], HEGRA [20], and EAS-TOP [21] experiments, and this work, 
as indicated in the legend. 
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